Influences of spin accumulation on the intrinsic spin Hall effect in two dimensional 
electron gases with Rashba spin-orbit coupling 
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In a two dimensional electron gas with Rashba spin-orbit coupling, the external electric field may 
cause a spin Hall current in the direction perpendicular to the electric field. This effect was called 
the intrinsic spin Hall effect. In this paper, we investigate the influences of spin accumulation on 
this intrinsic spin Hall effect. We show that due to the existence of boundaries in a real sample, 
the spin Hall current generated by the intrinsic spin Hall effect will cause spin accumulation near 
the edges of the sample, and in the presence of spin accumulation, the spin Hall conductivity will 
not have a universal value. The influences of spin accumulation on the intrinsic spin Hall effect in 
narrow strips of two dimensional electron gases with Rashba spin-orbit coupling are investigated in 
detail. 

PACS numbers: 72.10.-d, 72.15.Gd, 73.50.Jt 



I. INTRODUCTION 



Spintronics, which aims at the manipulation of the 
electron's spin degree of freedom in electronic devices, has 
become an emerging field of condensed matter physics 
for its potential application in information industryiS*^ 
Though originally many spintronic concepts involve fer- 
romagnetic metals since spins in this kind of sys- 
tem behave collectively and hence are easier to be 
controUed^i^iSi^ spintronics in semiconductors is more 
interesting because dope and heterojunction formation 
can be used to obtain some specific devices. It is an- 
ticipated that combining the advantages of semiconduc- 
tors with the concepts of spintronics will yield fascinat- 
ing new electronic devices and open the way to a new 
field of physics, i.e., semiconductor spintronics. How- 
ever, at present many great challenges still remain in 
this exciting quest. Among them, an issue that is fun- 
damentally important in semiconductor spintronics and 
has not been resolved is how to achieve efficient injec- 
tions of spins into non-magnetic semiconductors at room 
temperatureSiiSiiiiiSiiSi. Usage of ferromagnetic metals 
as sources of spin-injection is not practical because most 
of the spin-polarizations will be lost at the interface be- 
tween metal and semiconductor due to the large con- 
ductivity mismatch2ii£. Another possible approach is to 
use ferromagnetic semiconductors ( such as Gai_a;Mnj:As 
) instead of ferromagnetic metals as sources of spin- 
injection. In this approach, the problem of conductiv- 
ity mismatch does not exist and hence efficient injections 
of spins into non-magnetic semiconductors can truly be 
achievediiiiSii^. But for practical use at room tempera- 
ture, the Curie temperatures of ferromagnetic semicon- 
ductors are still too low. Thus, from both the exper- 
imental and theoretical points of view, more great ef- 
forts are still needed in order to achieve efficient injec- 
tions of spins in non-magnetic semiconductors at room 
temperature. Recently, a surprising effect was predicted 



theoretically that an electric field may cause a quantum 
spin Hall current in the direction perpendicular to the 
electric field in conventional hole-doped semiconductors 
( such as Si, Ge, and GaAs or in two dimensional elec- 
tron gases ( 2DEGs ) with Rashba spin-orbit couplingii^ 
This intrinsic spin Hall effect might offer a new approach 
for achieving efficient injections of spins in non-magnetic 
semiconductors and reveal a new avenue in the spintron- 
ics research. 

In this paper, we study the influences of spin accumu- 
lation on the intrinsic spin Hall effect in two dimensional 
electron gases with Rashba spin-orbit coupling. From 
the standpoint of spintronic applications, it is important 
to understand whether the spin Hall currents predicted 
in Refs. [14-15] are transport spin currents, i.e., whether 
they can be employed for transporting spins. An impor- 
tant feature of transport spin currents is that they will 
induce nonequilibrium spin accumulation at some spe- 
cific locations, for example, at the boundaries of a sample 
or at the interfaces between two different materialsSi^. 
On the other hand, if spin accumulation was caused in 
a sample due to the flow of a spin current, the spin 
current will also be changed significantly by the spin 
accumulation. ^'^'^ So in a strict theoretically treatment 
of the intrinsic spin Hall effect, the interplay between 
the spin Hall current and the spin accumulation must be 
taken into account. A detailed theoretical investigation 
of the influences of spin accumulation on the intrinsic 
spin Hall effect in hole-doped semiconductors was pre- 
sented in Ref.[16]. In the present paper, we will use a 
similar method as was applied in Ref.[16] to investigate 
the influences of spin accumulation on the intrinsic spin 
Hall effect in two dimensional electron gases with Rashba 
spin-orbit coupling. We will show that in contrast with 
what was found in Refs. [15, 17-18], in the presence of spin 
accumulation, the spin Hall conductivity in the intrinsic 
spin Hall effect in a Rashba two dimensional electron gas 
does not have a universal value, and in order to calculate 
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correctly the spin Hall current and the spin Hall conduc- 
tivity in a real sample with boundaries, the influences of 
spin accumulation need to be taken into account. The 
paper will be organized as follows: In Sec. II, we will first 
present a brief introduction to the intrinsic spin Hall ef- 
fect in a Rashba two dimensional electron gas in the ab- 
sence of spin accumulation and impurity scattering. In 
Sec. Ill, the influences of spin accumulation and impurity 
scattering will be taken into account. In Sec. IV, by use 
of the formulas derived in Sec. Ill, the influences of spin 
accumulation on the intrinsic spin Hall effect in a narrow 
strip of a Rashba two dimensional electron gas will be 
discussed in detail. 



II. INTRINSIC SPIN HALL EFFECT IN THE 
ABSENCE OF SPIN ACCUMULATION AND 
IMPURITY SCATTERING 

In this paper, we will use a slightly different method 
from what was applied in Refs. [15, 17-18] to discuss the 
intrinsic spin Hall effect in two dimensional electron 
gases with Rashba spin-orbit coupling. The merit of this 
method is that the influences of spin accumulation can 
be easily included. For clarity, in this section we first 
present a brief introduction to the intrinsic spin Hall ef- 
fect in the absence of spin accumulation and impurity 
scattering. We will show that our method will yield the 
same results as was obtained in Refs. [15, 17-18] if spin 
accumulation and impurity scattering are neglected. 

In the momentum representation, the single-particle 
Hamiltonian for a two dimensional electron gas with 
Rashba spin-orbit coupling readai^, 

Ho = — Q!(z X k) • cr 
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2m 



ak{cos9ay — sin^cra,). 



(1) 



where z is the unit vector along the normal of the two- 
dimensional plane, k = (fc cos 9, k sin 6) is the momentum 
of an electron, a —{<Jx, <^yi(^z) the Pauli matrices, and a 
the Rashba spin-orbit coupling constant. The Hamilto- 
nian (1) can be diagonalized exactly. The eigenenergies 
are given by 
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(A = ±l), 



(2) 



and the corresponding spinor eigenstates are defined by 
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(3) 



The expectation value of the spin of an electron which is 
in the spinor eigenstate |kA) will be given by 



where and Sy are the unit vectors along the x and y 

axes, respectively, and 5'^°],(k) and 5'^°^(k) are the x and 
y components of the spin, which are given by 
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5i°i(k) = |(kA|a,|kA) 



2k ' 



-A- 



2k 



(5) 
(6) 



Eqs.(4-6) show that the eigenstates of the Hamiltonian 
(1) are spin-polarized in the in-plane directions, and the 
spin-polarization directions depend on the momentum k. 

When an external electric field is applied in the x di- 
rection, the total Hamiltonian of the system will be given 
hy H = Hq + eE • X, where E =Exex is the external elec- 
tric field and — e the charge of the electron. The equation 
of motion for the electron's position and spin degrees of 
freedom under the action of the external electric field can 
be obtained directly from the Heisenberg equation, and 
the results read 



dt 
I. 



eE 

Kff X S, 



(7) 

(8) 



where he// = ^{kye^ — k^Gy) is an effective magnetic 
field felt by an electron with momentum k due to the 
Rashba spin-orbit coupling. The time variation of the 
spin of an electron which is initially in the spinor eigen- 
state |kA) ( i.e., S{t = 0) =|(kA|CT|kA) ) can be got by 
solving Eqs.(7-8) simultaneously. In this paper we will 
consider only the linear response of the transport prop- 
erty to the electric field. In the linear response regime, 
Eqs.(7-8) can be integrated analytically by use of the 
same method of Ref.[15], and one can find that under 
the action of the external electric field, the spin of an 
electron with momentum k and initially in the spinor 
eigenstate |kA) will become 



^A,.(k) 



ehky 



(9) 
(10) 

(11) 



where iS'A,i(k) is the zth component of the spin. Eqs.(9- 
11) show that an applied electric field in the x direction 
will cause the spin to tilt in the perpendicular direction 
by an amount proportional to ky. Due to this fact, the 
application of an external electric field in the x direction 
will induce a spin Hall current in the y direction with 
spin parallel to the z direction. The spin Hall current 
can be calculated by the following formula 



5A,.(k)(;ifc,/m)/A(k), (12) 



Sf(k) 



-(kA|a|kA) 

^S(k)e. + 5S(k)e, 



(4) 



where /A(k) is the probability distribution function of 
conduction electrons. If spin accumulation and electron- 
impurity scattering can be neglected, f\ (k) can be given 
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simply by the Fermi-Dirac equilibrium distribution func- 
tion, i.e., Mk) = /O(ekA) = l/{exp[/3(ekA - ep)] + 1}, 
where /3 = l/fc^T and ep is the Fermi energy. Then 
one can find that in the usual case where both spin-orbit 
split bands are occupied, the spin Hall current Jy' and 
the spin Hall conductivity (Ts will be given hy^ 



-'y 



f^sEx 



(13) 



Eq.(13) is the central result of Ref.[15], which was also 
obtained by several other groups with different theoret- 
ical approachesfii Eq.(13) shows that as long as both 
spin-orbit split bands are occupied, the spin Hall conduc- 
tivity will have a universal value, independent of both 
the Rashba spin-orbit coupling strength and of the den- 
sity of conduction electrons. Though Eq.(13) was ob- 
tained in the clean limit, recent numerical simulation 
shows that it still holds in the presence of weak disor- 
der, providing that the sample size exceeds the localiza- 
tion length.^* It is also important to note that unlike the 
similar effect conceived by Hirscb22., which is caused by 
spin-orbit dependent anisotropic scattering from impuri- 
ties and will vanish in the weak scattering limit f22i2iiSS 
the spin Hall effect described by Eq.(13) has a quantum 
nature and is purely intrinsic, i.e., it does not rely on 
any anisotropic scattering from impurities. Of course, 
it should be pointed out that though the mechanism of 
the intrinsic spin Hall effect described above does not 
involve impurity scattering, it does not mean that impu- 
rity scattering have no significant influences on the effect. 
The reason is that in a real sample, due to the existence 
of boundaries, nonequilibrium spin accumulation will be 
caused inevitably near the edges of the sample when the 
spin Hall current circulates in it, and in the presence of 
nonequilibrium spin accumulation, spin diffusion will be 
induced by electron-impurity scattering and, hence, the 
spin Hall current may also be changed significantly from 
what was given by Eq.(13). Thus in order to calculate 
correctly the spin Hall current and the spin Hall conduc- 
tivity in a real sample with boundaries, the influences of 
spin accumulation and electron-impurity scattering need 
to be considered. 



III. INFLUENCES OF SPIN ACCUMULATION 
AND IMPURITY SCATTERING 

In the presence of spin accumulation and electron- 
impurity scattering, the distribution function of conduc- 
tion electrons can no longer be given simply by the Fermi- 
Dirac equilibrium distribution function but should be de- 
rived strictly by solving the Boltzmann transport equa- 
tion. In a steady ( but nonequilibrium ) state, the Boltz- 
mann equation reads 



VA(k).V/A(r,k)-eEe,fVA(k) 



'9/A(r,k) 



[(^)i^.T^^ + (^)i^.T^^], (14) 



where VA(k) = |;VkekA is the velocity of conduction 
electrons, Eext — E^e^is the external electric field ap- 
plied in the x direction, and /^(r, k) is the distribution 

function. The collision term [dfx/dt) 



(A^A ) 
coll. 



describes the 

changes of the distribution function due to the intra-band 
( A = A ) and/or inter-band (A = A ) electron-impurity 
scattering, which is given by 



\ 5t Icon 



X [/A(r,k)-/v(r,k')], 



(k,k')<5(eA(k)-6v(k')) 
(15) 



where u)^ y (k, k') is the rates of an electron to be scat- 
tered from the state |kA) into the state |k'A ) by impurity 
scattering. 

The Boltzmann equation (14) can be solved by the 
relaxation time approximation method. Within the re- 
laxation time approximation and in the linear response 
regime, the system can be considered only slightly de- 
viated from the equilibrium state, thus the total dis- 
tribution />(r,k) can be expressed as the sum of the 
equilibrium distribution function /"(eka) and the non- 
equilibrium ones as the following 



/A(r,k) =/«(ekA)-eM\r) 



dfie^^x) 



9ekA 



erA(k)E^(r).VA(k) 



dekx 



(16) 



Here the second term denotes the change of the distri- 
bution function due to the occurrence of nonequilibrium 
spin accumulation in the sample, with — e/Lt^(r) ( A = ± 
) denoting the band- and position-dependent shifts of 
the Fermi level in the nonequilibrium state, which char- 
acterize the imbalance of the filling of conduction elec- 
trons in the two spin-orbit split bands in the presence 
of nonequilibrium spin accumulation. The third term in 
Eq.(16) denotes the changes of the distribution function 
due to the movement of conduction electrons under the 
action of an effective electric field E'*'(r), with r;^(k) de- 
noting the total relaxation time of conduction electrons 
with momentum k, which is determined by the electron- 
impurity scattering. Because the occurrence of nonequi- 
librium spin accumulation will cause spin diffusion in the 
sample, the effective electric field E'^(r) should be given 
by E'^(r) = E^xt — V/i'^(r), i.e., in addition to the exter- 
nal electric field Eext, conduction electrons will also feel 
an effective field given by the gradients of the band- and 
position-dependent shifts of the Fermi level&LS Substi- 
tute Ea. l(TC|l into Eqs. (14-15), the left-hand side of the 
Boltzmann equation will become 



Lhs 



-e^4^{E^r).VA(k) 
-n(k)VA(k).v[E^r).VA(k)]}, (17) 



and the right-hand side of the Boltzmann equation be- 
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come 



nns = — e — {rA(kj[ 



(k) 



(18) 



where Tj^(k) and T|^(k) are the intra-band and inter- 
band transition relaxation times, respectively, which are 
defined by 

rri^) = [J ^c.A,A(k,k')<5(e.(k)-6,(k'))]-i (19) 

riUk) = [ J ^c..,_.(k,k')-5(6.(k) - 6_,(k'))]-^ 

(20) 

By integrating both the left-hand side and the right-hand 
side of the Boltzmann equation, one can find that the 
total relaxation time t\ of conduction electrons should 
be given by 



TA(k) 



1 



1 



(k) 



-Ti 



(k)^ 



(21) 



and the band- and position-dependent shifts of the Fermi 
level satisfy the following equation, 



,2,A.^ M^(i-)"M"^(r) 



vV"(r) 



(22) 



where Dx = [(V^)^t;-t-],^]^/^, with denoting the 
band-dependent Fermi velocity and Tp^ the inter-band- 
transition relaxation time and Tp the total relaxation 
time of conduction electrons at the Fermi level, respec- 
tively. ( For simplicity, we assume that Tp and Tp^ are 
band-independent. ). From Eq.(22), one can see that 
the relative shifts of the Fermi level in the two spin-orbit 
split bands, given by /i^(r) — A' (i")i satisfy the following 
diffusion equation, 

V^[,-(r)-,-(r)]^^^^i£b^, (23) 



where D is the spin-diffusion length, defined by 



D=\— + —1-1/2 



(24) 



In addition to Eqs. (22-23), the band-dependent shifts of 
the Fermi level should also satisfy the charge neutral- 
ity condition, which requires that the net changes of the 
charge density due to the band-dependent shifts of the 
Fermi level, given by Sp{r) = c^^a/ (f;;y2 [/a(i", k) - 

/^(ekA)]j should be zero. This requirement arises from 
the fact that according to the symmetry of the Hamilto- 
nian (1), in the direction perpendicular to the external 
electric field, no charge Hall current will be generated. 



so the occurrence of spin accumulation due to the flow 
of the spin Hall current does not result in charge accu- 
mulation. Due to this requirement, one can show that in 
addition to Eqs. (22-23), the band-dependent shifts of the 
Fermi level should also satisfy the following equation, 



J2 kxFfi\r) = 0, 



(25) 



A=± 



where k^p is the band-dependent wave number at the 
Fermi level. After the band-dependent shifts of the Fermi 
level are determined, the spin Hall current can be ob- 
tained by inserting Eq.(ll) and Eq.(16) into Eq.(12), 
then one can find that in the usual case where both spin- 
orbit split bands are occupied^^, the spin Hall current 
and the spin Hall conductivity will be given by 



(26) 
■•(27) 



Stt l6TTa[{ma/tfy + 2mep/h^]^/^■ 
Eqs.{26-27) show that the spin accumulation may have 
some significant infiuences on the intrinsic spin Hall effect 
in a Rashba two dimensional electron gas. Firstly, in the 
presence of spin accumulation, the spin Hall conductivity 
will be a position — dependent quantity and do not have 
a universal value, i.e., it will depend on the Rashba spin- 
orbit coupling constant and on the density of conduction 
electrons. This is very different from what was shown in 
Eq.(13). Secondly, in the presence of spin accumulation, 
the spin Hall current may be decreased substantially from 
the corresponding value obtained in the absence of spin 
accumulation, and the decrease will be determined by 
/i+(r) — /i~(r), i.e., proportional to the relative shifts of 
the Fermi level in the two spin-orbit split bands. 



IV. RESULTS AND DISCUSSIONS 

Eqs. (23,25) and (26-27) constitute a set of self- 
consistent equations, from which both the spin Hall cur- 
rent and the spin accumulation can be obtained. In 
this section, we apply these formulas to discuss the in- 
trinsic spin Hall in a narrow strip of a two dimensional 
electron gas with Rashba spin-orbit coupling. Narrow 
strips are the usual geometry applied in the experimen- 
tal measurement of the Hall effect, including the spin 
Hall effect In the following we will assume that 

the longitudinal direction of the strip is along the x axis 
and the transverse direction along the y axis and the nor- 
mal of the 2D plane along the z axis, respectively, and 
an external electric field is applied in the longitudinal 
direction of the strip. According to Eqs.(26)-(27), in or- 
der to calculate the spin Hall current Jy' caused by the 
longitudinal external electric field , one must first find 
out the band-dependent shifts fJ.~^{r) and of the 

Fermi level. For simplicity, we assume that the length L 
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of the strip is much larger than its width w so that spin 
diffusion in the longitudinal direction of the strip can be 
neglected. In such case, only transverse spin accumula- 
tion need to be considered, and /i'*'(r), Jy'{r), and crs(r) 
will all depend only on the y coordinates. From Eq. (23) , 
— ^J.^{y) can be expressed as 



(28) 



where A and B are two constant coefBcients that need 
to be determined by appropriate boundary conditions. 
In this paper, we will consider the transverse-open- 
circuit boundary condition. In the transverse-open- 
circuit boundary condition, the spin Hall current at the 
two boundaries of the sample, which are assumed to be 
located at y = i^, should be zero, i.e. 



±-) 



0. 



(29) 



Substituting Eq.lI2Hl) into Eqs. (26-27) and by use of the 
above boundary condition, the coefficients A and B can 
be determined. Then the band- and position-dependent 
shifts of the Fermi level in the strip can be obtained, and 
we get 



Ai"(2/) 



Mocosh(y/i^) D'' 
2cosh(w;/2L>) ^ Dl 



Mocosh{y/D)^D^_ 



2cosh(w/2D) 'Dl 
where Mq and Dq are defined by 



Mo 



Dn = 



— ^(ma/fi,2)2 + 2meF/h'^, 
e 

hi 



1-1/2 



(30) 
(31) 

(32) 
(33) 



Here D\ ( A = ± ) has been defined in Eq.(22). In ob- 
taining Eqs. (30-31), the charge neutrality condition (25) 
was also used. After ii^{y) is determined, according to 
Eqs. (26-27), the spin Hall current and the spin Hall con- 
ductivity will also be obtained, and the results read 



e cosh{y/D) 
Stt cosh(?i;/2L') 



(34) 
(35) 



Eqs. (34-35) show that, due to the influences of spin accu- 
mulation, the spatial distribution of the spin Hall current 
in a sample will be highly inhomogeneous and the spin 
Hall conductivity is sensitively position-dependent. The 
spin Hall current and the spin Hall conductivity and their 
spatial distributions will also have sensitive dependences 
on the spin diffusion length D and the sample width w. 
This was shown in Fig.l, where we have plotted the trans- 
verse spatial distributions of the spin Hall currents in 
three distinct cases with different ratios of D/w. From 
Fig.l, one can see that if w D, the spin Hall current 



will be negligibly small in the sample, i.e., Jy''{y) ~ 
everywhere. On the other hand, if w 3> -D, the spin 
Hall conductivity will be approximately a constant at 
|y| < w/2, i.e., asiy) = J^'{y)/E, =^ ^ at \y\ < w/2, 
which is independent of both the Rashba spin-orbit cou- 
pling strength and of the density of conduction electrons. 
But <Ts{y) will decrease to zero as y ^ ±u;/2. 

The spin accumulation caused by the longitudinal elec- 
tric field Ex can be calculated through the following for- 
mula 



E 



SA(k)/A(r,k), 



(36) 



where SA(k) has been given in Eqs. (9-11). By inserting 
Eq.(16) and Eqs.(9-ll) into Eq.(36) and with the help of 
Eqs. (30-31), one can find that both the y component of 
(S) ( the in-plane spin accumulation ) and the z compo- 
nent of (S) ( the perpendicular spin accumulation ) are 



non-zero, 



{Sy) 



emaTpEx 
eEx , rriTF 



,1/2(1 



a^m smh{y/D) 



h^ep cosh{wl2D) 



(37) 
.(38) 



Eqs. (37-38) show that both the in-planc and the perpen- 
dicular spin accumulation are proportional to the longi- 
tudinal electric field E^, but there are some significant 
differences between them. The in-plane spin accumu- 
lation is homogeneously distributed in the sample and 
independent of both the spin diffusion length D and of 
the sample width w. However, the spatial distribution of 
the perpendicular spin accumulation is highly inhomoge- 
neous and its magnitude depends sensitively on the spin 
diffusion length D and on the sample width w. These dif- 
ferences arise from the fact that the in-plane and the per- 
pendicular spin accumulation are caused by very different 
mechanism. In fact, it was known long time ago that in 
a two dimensional electron gas with Rashba spin-orbit 
coupling, an applied in-plane electric field will induce a 
homogeneous in-plane spin accumulation polarized in the 
direction perpendicular to the electric fieldS^iS^, but the 
in-plane spin accumulation has nothing to do with the in- 
trinsic spin Hall effect. From the theoretical viewpoints, 
the in-plane spin accumulation is caused by the combined 
action of the spin-orbit coupling, absence of inversion 
symmetry, and the time-reversal symmetry-breaking in 
the electric fieldSiS^. Since the in-plane spin accumula- 
tion has been investigated in detail in previous literatures 
and it has no relation with the intrinsic spin Hall effect, 
we will not discuss it again in the present paper. Unlike 
the in-plane spin accumulation, the perpendicular spin 
accumulation given by Eq.(38) is caused by the intrinsic 
spin Hall effect, so its spatial distribution is highly inho- 
mogeneous and its magnitude depends sensitively on the 
spin diffusion length D and on the sample width w. This 
was illustrated clearly in Fig. 2, where we have plotted 
the transverse spatial distributions of the perpendicular 
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spin accumulation in three distinct cases with different 
ratios oi D/w. From Fig. 2 and Eq.(38), one can see that 
the perpendicular spin accumulation is maximum at the 
edges of the sample, and the perpendicular spin accumu- 
lation at the edges of the sample will increase with the 
increase of the width of the sample. When the width 
w of the sample is much larger than the spin diffusion 
length D, the perpendicular spin accumulation at the 
edges of the sample will approach a maximum value of 
^i^TT^y^'^i^+WT^)^ which is independent of the sam- 
ple width w. This will be a merit for the experimental 
measurement of the effect. In order to get a quantita- 
tive estimation of the perpendicular spin accumulation, 
let us consider some actual experimental parameters. In 
current 2DEG high quality samplesi^ii^iS^, the typical 
carrier concentrations range from 5 x 10^^ to 10^^ cm~^, 
the strength of the Rashba spin-orbit coupling is on the 
order of 1 x 10^^"'^ 5 x 10~^^eVm, the effective mass 
of conduction electrons is about 0.05me, the relaxation 
time is typically Ips, the spin diffusion length is about 
Ijim, and the Fermi energy e_F is about 20 ~ 50meV. If 
one consider a sample with the width w = 10/xm ( much 
larger than the spin diffusion length ) and the external 
field eEx = lOKeV/m, then from Eq.(38) one can es- 
timate that the perpendicular spin accumulation at the 
edges of the sample can be as larger as 10~^^J • sjm} . 
This magnitude should be large enough to be detected 
exp eriment ally. 




FIG. 1: Illustration of the transverse spatial distributions of 
the spin Hall currents in three distinct cases with different 
ratios of D/w. The parameters used are: the sample width 
w = lO^m; the spin diffusion length D — Ifim (the solid line), 
10/im (the dashed line), and 100pm (the dotted hne). 

In conclusion, in this paper we have investigated the in- 
fluences of spin accumulation on the intrinsic spin Hall ef- 
fect in two dimensional electron gases with Rashba spin- 
orbit coupling. We have presented a detailed theoretical 
analysis on the interplay between the spin Hall current 
and spin accumulation in the intrinsic spin Hall effect in 
a Rashba two dimensional electron gas. We have shown 
that in the presence of spin accumulation, the spin Hall 
conductivity will not have a universal value. The spin 
Hall current and spin accumulation in narrow strips of 




FIG. 2: The transverse spatial distributions of the perpendic- 
ular spin accumulation in three distinct cases with different 
ratios of D/w. ( The parameters used are: the sample width 
w — 10 /im; the spin diffusion length D — 1pm (the solid 
line), 10/xm (the dashed line), and 100/im (the dotted line); 
the external field eE^ ~ lOKeV/m; the Rashba spin-orbit 
coupling constant q = 1 x 10~^^eVm; the effective mass m = 
0.05me; the relaxation time is Ips; and the Fermi energy ep = 
20meV.) 



two dimensional electron gases with Rashba spin-orbit 
coupling was calculated explicitly. The results show that 
in order to calculate correctly the spin Hall current and 
the spin Hall conductivity in a real sample with bound- 
aries, the influences of spin accumulation need to be 
taken into account. Recently, E. I. Rashba pointed out 
that the Hamiltonian (1) implies that there exist non- 
vanishing dissipationless spin currents even in the ther- 
modynamic equilibrium state ( i.e., in the absence of the 
external electric field ) ^ These background spin currents 
are not associated with real spin transports but spurious 
effects caused by the lacking of the time-reversal sym- 
metry implied in the Hamiltonian (1). Due to this fact, 
a procedure for eliminating the spurious effects of these 
background spin currents should be devised in calculat- 
ing transport spin currents if the background currents 
contribute to the calculation. But for the intrinsic spin 
Hall effect discussed in the present paper, the background 
spin currents do not contribute to the calculation of the 
spin Hall current due to the following reasons. First, the 
spin Hall current is polarized in the direction perpendic- 
ular to the 2D plane, while the background spin currents 
are polarized in the 2D plane. Second, the spin Hall cur- 
rent is a dynamic response of the spins to the external 
electric field and will vanish in the absence of the electric 
field, but the background spin currents are independent 
of the electric field. Due to these reasons, the background 
spin currents do not present in the calculation of the spin 
Hall current and hence don't need to be considered in the 
present paper. 
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